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1.  Introduction 

Populations  of  cells  that  make  up  organ  tissue  grow  and  contract  in  a  manner 
that  is  roughly  analogous  to  the  fluctuations  of  other  natural  populations.  Since 
organs  are  of  bounded  size  their  growth  is  not  entirely  uncontrolled  and 
exponential  overall,  or  otherwise  a  human  or  rat  liver  would  either  assume  a 
totally  outlandishly  large  size,  or  else  shrink  to  zero.  Neither  such  alternative  is 
seen  in  nature,  although  organ  sizes  do  vary  among  otherwise  comparable 
members  of  the  same  species.  And  organs  of  mature  hosts  can  change  in  size  as  a 
result  of  disease,  toxic  insult,  or  an  operation  such  as  partial  hepatectomy,  from 
which  a  liver  can  recover  again  to  normal  size  and  function. 


2 


A  traditional  approach  to  modeling  population  (=  organ)  size  restriction  to  an 
observed  "normal"  level  is  to  postulate  a  carrying  capacity :  effectively  a  limit  on 
the  physical  region  that  can  be  occupied  by  the  population.  This  formulation 
apparently  goes  back  to  Verhulst  (1836);  see  Murray  (1989)  for  recent  discussion. 
In  the  organ  situation  this  might  correspond  to  a  space  of  approximately 
pre-ordained  dimension  that,  say,  liver  cells  in  liver  tissue  cannot  exceed  in  the 
body  of  a  mature  human  male.  The  space  can  be  taken  as  given,  introduced  into 
other  models  as  a  parameter,  and  in  particular  cases  estimated  from  data.  It 
would  be  the  maximum  size  of  the  liver  compartment  in  a  PB-PK  model,  for 
example. 

There  is  another  alternative  to  the  above  approach  that  depends  upon 
recognition  and  measurement  of  the  presence  of  various  biological  agents  called 
growth  factors.  There  are  a  number  of  such  factors  that  both  stimulate  (positively) 
and  inhibit  (negatively)  cell  population  growth.  Growth  factors  are  discussed  by 
Alberts  et  al.  (1994).  Aaronson  (1991)  provides  an  overview  of  growth  factors  in 
cancer;  see  also  Rubin,  Bottaro,  and  Aaronson  (1993).  The  purpose  of  this  note  is 
to  provide  a  very  simple  model  for  a  cell  population  that  grows  under  the  control 
of  positive  and  negative  growth  factors.  It  will  be  seen  that  such  a  model  can 
result  in  logistic  growth  without  the  necessity  of  postulating  a  physical  carrying 
capacity.  An  effective  carrying  capacity  appears  as  a  result  of  presumed  growth 
factor  interaction  with  cells. 

2.  Model  for  a  Cell  Population  Under  Growth  Factor  Control 

Suppose  a  population  of  cells  is  of  size  C(f)  at  time  t.  Its  individual  cell  growth 
or  birth  rate  is  Ao,  and  its  death  rate  is  po,  so  its  net  growth  rate,  Ao  -  [1q,  governs 
the  manner  and  rate  of  growth.  Starting  with  C(0)  members,  and  left  alone,  the 
population  would  grow  roughly  like  C(£)  ~  C(0)e^~^1 ,  which  means  either  to 
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a  large  size  (Ao  -  IM)  >0),  or  to  zero  (Ao-/io<0)-  Clearly  such  unrestricted 
behavior  is  inappropriate  for  describing  a  population  of  cells  that  constitutes  an 
entire  organ,  although  essentially  such  a  model  has  been  used  to  describe  growth 
of  tumors  within  an  organ;  see  Tan  (1991)  for  an  overview;  in  particular  the  work 
of  Moolgavkar  and  co-authors,  cited  in  Tan  (1991),  is  relevant. 

Now  introduce  a  quantity  ait)  of  a  positive  growth  factor  into  the  vicinity  of  the 
cell  population.  The  amount  present,  ait),  changes  cell  birth  rate  to  Ao  +  X\  a(t), 
where  we  take  Ai  >  0.  Also  introduce  a  quantity  Pit)  of  negative  growth  factor;  it 
changes  cell  death  rate  to  po  +  Pipit),  p\  >  0.  Then  changing  levels  of  ait)  and  pit) 
can  certainly  alter  the  properties  of  the  cell  population,  from  net  growth  to  net 
decline,  depending  upon  values  of  a(t)  and  pit). 

Assume  that  the  productions  of  both  ait)  and  Pit)  are  regulated  by  cell 
activity  in  such  a  way  that 

Yt  =PaC(t)-5aa(t )  (2.1) 

and 

yf  =  ppC(t)-5pm-  (2.2) 

That  is,  both  are  stimulated  to  increase  by  the  number  of  cells  present,  and  to 
diminish  in  proportion  to  their  own  concentration,  possibly  being  removed  from 
the  cell  site  (organ)  by  blood  flow  or  metabolism  or  other  biological  processes.  Of 
course  the  above  equations  are  prime  candidates  for  replacement  by  others  that 
more  accurately  depict  the  true  interactions. 

In  the  presence  of  a  and  /^-factors  the  cells  in  the  organ  grow  and  decline 
according  to 

=  [A0  +  A1a(f)]C(f)  -  [^o  +  (2.3) 
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So  (2.1),  (2.2),  (2.3)  form  a  system  of  three  non-linear  differential  equations.  No 
explicit  solution  seems  immediately  available,  unless  we  make  the  quasi-static  or 
quasi-steady-state  assumption  (QSSA);  see  Strogatz  (1994)  for  its  invocation  so  as 
to  solve  a  non-linear  dynamics  problem  along  with  some  historical  references, 
and  Segel  and  Slemrod  (1989)  for  a  careful  discussion  of  this  approximation. 
Namely  assume  that  a(t )  and  (3(t )  are  able  to  adapt  very  quickly  to  any  current 
value  of  C(f)  to  always  reach  a  "temporary  steady  state": 

—  0  =  Pa^-(0  —  3aa(t)  (2.4) 

^  =  0  =  p/JC(t)-S/!/)(t).  (2.5) 

Adopt  the  approximation  as  true,  so  solve  (2.4)  and  (2.5)  for  a(t)  and  J3(t): 

*(*HPa/W  t)  ( 2.6 ) 

and 

Pit)  =  {pp/Sp)c(t).  (2.7) 

Let  us  call  (pa/  $a)  and  (pp/  Sp)  the  prevalences  of  the  a  and  (5  factors  respectively. 
Insert  these  into  (2.3)  and  for  convenience,  put  X{  =  ^i(pa/<5a),  p{  =  P\[pp /^p)> 

to  obtain 

=  [Ao  +  -  [ft,  +  tfC(( )]c(f ) 

=  (4-rt))C(()-(M-Al')C2W  (2.8) 

=  (A0-Wl)C(i{l-p^C(f). 

“Po 

This  conforms  exactly  to  the  original  logistic  equation  if  the  ordinary  net  growth 
rate,  A  =  -  po,  is  positive,  as  is  the  effective  carrying  capacity 


Under  the  above  conditions  and  starting  from  C(0)  >  0,  the  population  attains  the 
long-run  steady-state  value 

CH  =  K - V?°  .  (2.10) 

^1  [Pp  /dp )  M  (Pa  /d a ) 

The  above  version  of  carrying  capacity  makes  intuitive  sense  in  that 

(a)  it  increases  with  net  population  growth  rate,  Xo  -  po; 

(b)  it  decreases  with  increased  prevalence  of  negative  growth  factor,  {pp/ dp), 
and  with  decreased  prevalence  of  positive  growth  factor,  ip  a/  da)', 

(c)  the  inhibition  effect  of  negative  growth  factor,  p\  =  p\{pp/ dp),  must 
exceed  the  stimulative  effect  of  the  positive  growth  factor, 

X{  =  X]  {pa/  da). 

If  any  of  the  above  conditions  are  violated  the  population  development 
becomes  radically  different,  but  can  also  be  interesting. 

The  time-dependent  population  size  is  seen  to  be  of  the  familiar  logistic 
growth  form 


C(t)  = 


KC(0)el 


K-C(0)  +  C(0)eAt 

with  K  as  in  (2.10),  A  =  Xq  -  po  >  0,  and  0  <  C(0)  <  K. 


(2.11) 


Note  that  the  formula  has  biological  meaning  even  if  C(0)  >  K,  and  also  if 
K<0:  suppose  that  A  =  Xq  -  po  >  0  but  X[  >  p{;  then  write  K'  =  -K>  0  to  get 


^=AC(0[i+c(t)/r], 


(2.12) 


the  solution  to  which  is 


(2.13) 


K'C(  0)g 


Af 


K'  +  C(0)  -  C(0)e 


Af 


if  t<—  In 
A 


K'  ] 

,  1 

f  K'  1 

C(0)J 

;  it  explodes  when  t  =  —ln 

C(0)J 

;  this  might  plausibly 


model  an  especially  malignant  tumor  growth.  Finally  when  A  =  Aq  -  fio  <  0  and 
K'  =  -K  >  0  we  simply  get  (2.11)  with  -A  replacing  A,  once  again  a  logistic  model, 
but  now  one  that  decreases  as  t  increases. 


3.  Stochastic  Models 

The  above  model  can  be  "made  stochastic"  in  various  ways.  One  is  to  re-state 
the  growth  factor  and  cell-growth  equations  as  a  system  of  three  non-linear 
Ito-type  stochastic  differential  equations.  Analytical  solutions  are  not  likely  to  be 
available,  but  some  asymptotics  might  well  produce  explicit  results. 

Another  approach  is  to  assume  that  a  stochastic  version  of  C(t),  namely  C(f),  is 
a  birth-and-death  process  with  transition  rates  copied  from  the  right-hand  side  of 
(2.8).  For  example,  C(f)  evolves  over  state  space  (0, 1, 2, ...)  according  to 

P{c(t  +  dt)  =  C(t)  +  l|C(f)  =  i }  =  (Aq-  fio)idt  +  o(dt )  (3.1,a) 

•  p{C(t  +  dt)  =  C(f)  - l|C(f)  =  i]  =  (n[  -  X{)i2dt  +  o{dt);  (3.1,b) 

these  holding  when  A  =  Ao  -  no  >  0,  [i{  -  X{  >  0;  otherwise  modification  is  needed. 
Although  £[C(f)|C(0)  =  C(0)j  ^  C(t )  because  of  the  non-linearity  in  the  generator, 

(3.1),  it  can  be  of  interest  to  study  the  above  stochastic  version's  transient 
properties,  such  as  first-passage  times  from  low  states  (population  size)  to  high, 
or  the  reverse,  e.g.  to  C(f#)  =  0  when  the  population  dies  out. 

The  described  approach  essentially  minimizes  attention  to  the  stochastics  of  the 
growth  factors  and  ignores  non-linearity,  hence  is  a  prime  candidate  for  an 
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upgraded  treatment.  Nevertheless  it  is  appealing  for  its  simplicity  and  easy 
availability,  and  is  offered  as  an  interim  approach. 

4.  Summary 

It  is  shown  that  classical  logistic  growth  can  be  induced  in  a  non-traditional 
manner  by  hypothesized  action  of  growth  factors,  rather  than  by  action  of  a 
physical  carrying  capacity  (although  the  latter  may  operate  as  well).  Modification 
of  the  effective  carrying  capacity  to  be  negative  has  biological  interpretability. 
The  resulting  models  may  perhaps  find  a  use  in  cell  proliferation  and  cancer 
modeling. 
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29.  Dr.  Andrew  Gelman . 1 

Statistics  Dept. 

University  of  California 
Berkeley,  CA  94720 


30.  Dr.  NeilGerr . 1 

Office  of  Naval  Research 
Arlington,  V A  22217 
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31.  Prof.  Peter  Glynn . 

Dept,  of  Operations  Research 

Stanford  University 
Stanford,  CA  94305 

32.  Prof.  Linda  V.  Green . 

Graduate  School  of  Business 
Columbia  University 

New  York,  NY  10027 

33.  Prof.  J.  Michael  Harrison . 

Graduate  School  of  Business 

Stanford  University 
Stanford,  CA  94305-5015 

34.  Dr.  D.  C.  Hoaglin . 

Department  of  Statistics 

Harvard  University 
1  Oxford  Street 
Cambridge,  MA  02138 

35.  Dr.  David  G.  Hoel . 

Professor  of  Biometry  and  Epidemiology 
Medical  University  of  South  Carolina 

171  Ashley  A ve. 

Charleston,  SC  29425-0002 

36.  Prof.  D.  L.  Iglehart . 

Dept,  of  Operations  Research 

Stanford  University 
Stanford,  CA  94305-5015 

37.  Institute  for  Defense  Analysis . 

1800  North  Beauregard 

Alexandria,  VA  22311 

38.  Dr.  Robert  C.  Jackson . 

Vice  President,  Research  and  Development 
Agouron  Pharmaceuticals,  Inc. 

3565  General  Atomics  Court 
San  Diego,  CA  92121-1121 

39.  Prof.  J.  B.  Kadane . 

Dept,  of  Statistics 

Camegie-Mellon  University 
Pittsburgh,  PA  15213 

40.  Dr.  F.  P.  Kelly . 

Statistics  Laboratory 

16  Mill  Lane 
Cambridge,  ENGLAND 
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41.  Dr.  Jon  Kettenring . 1 

Bellcore 

445  South  Street 

Morris  Township,  NJ  07962-1910 

42.  Dr.  Ralph  Kodell . 1 

Chief,  Biometry  Branch 

Biometry  and  Risk  Assessment 


National  Center  for  Toxicology  Research 
3900  NCTR  Drive 
Jefferson,  AR  72079 

43.  Mr.  Koh  Peng  Kong . 1 

OA  Branch,  DSO 

Ministry  of  Defense 
Blk  29  Middlesex  Road 
SINGAPORE  1024 

44.  Prof.  Guy  Latouche . 1 

University  Libre  Bruxelles 

C.P.  212,  Blvd.  De  Triomphe 
Bruxelles  B-1050,  BELGIUM 

45.  Dr.  A.  J.  Lawrance . ^ . 1 

Dept,  of  Mathematics 

University  of  Birmingham 
P.O.  Box  363 

Birmingham  B15  2TT,  ENGLAND 

46.  Prof.  M.  Leadbetter . 1 

Department  of  Statistics 

University  of  North  Carolina 
Chapel  Hill,  NC  27514 

47.  Prof.  J.  Lehoczky . 1 

Department  of  Statistics 

Camegie-Mellon  University 
Pittsburgh,  PA  15213 

48.  Dr.  Georg  Luebeck . 1 

Fred  Hutchinson  Cancer  Research  Center 

1124  Columbia 
MP-665 

Seattle,  WA  98014 

49.  Dr.  Colin  Mallows . 1 

AT&T  Bell  Telephone  Laboratories 

600  Mountain  Avenue 
Murray  Hill,  NJ  07974 
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1 


50.  Prof.  R.  Douglas  Martin . . . 1 

Department  of  Statistics,  GN-22 

University  of  Washington 
Seattle,  WA  98195 

51 .  Dr.  Sati  Mazumdar . .  \ 

Biostatistics  Dept. 

University  of  Pittsburgh 
Graduate  School  of  Public  Health 
Pittsburgh,  PA  15261 

52.  Dr.  James  McKenna . . . .  i 

Bell  Communications  Research 

445  South  Street 
Morristown,  NJ  07960-1910 

53.  Dr.  Ramit  Mehr-Grossman . . . 1 

Theoretical  Biology  and  Biophysics 

Theoretical  Division 
Mail  Stop  K710 

Los  Alamos  National  Laboratory 
Los  Alamos,  NM  87545 

54.  Prof.  Carl  N.  Morris . . . . . . . . . .  \ 

Statistics  Department 

Harvard  University 
1  Oxford  Street 
Cambridge,  MA  02138 

55.  Dr.  John  A.  Morrison . . . . . \ 

AT&T  Bell  Telephone  Laboratories 

600  Mountain  Avenue 
Murray  Hill,  NJ  07974 

56.  Prof.  F.  W.  Mosteller . . . \ 

Department  of  Statistics 

Harvard  University 
1  Oxford  Street 
Cambridge,  MA  02138 

57.  Dr.  John  Orav . \ 

Biostatistics  Department 

Harvard  School  of  Public  Health 
677  Huntington  Avenue 
Boston,  MA  02115 

58.  Dr.  Alan  Perelson . \ 

Theoretical  Biology  and  Biophysics 

Theoretical  Division 
Mail  Stop  K710 

Los  Alamos  National  Laboratory 
Los  Alamos,  NM  87545 


1 


59.  Dr.  Jim  Petty . 

National  Biological  Survey 
4200  New  Haven  Road 
Columbia,  MO  65201 

60.  Dr.  Lorenz  Rhomberg . 1 

Harvard  Center  for  Risk  Analysis 

Harvard  University 
Cambridge,  MA  02138 


61 .  Dr.  Rhonda  Righter . 1 

Dept,  of  Decision  &  Info.  Sciences 

Santa  Clara  University 
Santa  Clara,  CA  951 18 

62.  Dr.  John  E.  Rolph . 1 


Information  and  Operations  Management 
Univ.  of  Southern  California 
School  of  Business  Administration 
Los  Angeles,  CA  90089-1421 

63.  Prof.  M.  Rosenblatt . 1 

Department  of  Mathematics 

University  of  California,  San  Diego 
La  Jolla,  CA  92093 

64.  Prof.  Frank  Samaniego . 1 

Statistics  Department 

University  of  California 
Davis,  CA  95616 

65.  Prof.  G.  A.  F.  Seber . 1 

Dept,  of  Statistics 

Univ.  of  Auckland 
Private  Bag  92019 
Auckland,  NEW  ZEALAND 

66.  Prof.  G.  Shantikumar . 1 

The  Management  Science  Group 

School  of  Business  Administration 
University  of  California 
Berkeley,  CA  94720 


67.  Prof.  N.  D.  Singpurwalla . 1 

George  Washington  University 

Washington,  DC  20052 

68.  Prof.  H.  Solomon . 1 

Department  of  Statistics 

Sequoia  Hall 
Stanford  University 
Stanford,  CA  94305 
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69.  Dr.  Andrew  Solow . 

1 

Woods  Hole  Oceanographic  Institute 

Woods  Hole,  MA  02543 

70.  Prof.  W.  Stuetzle . . . 

1 

Department  of  Statistics 

University  of  Washington 

Seattle,  WA  98195 

71.  Prof.  J.  R.  Thompson . 

1 

Dept,  of  Mathematical  Science 

Rice  University 

Houston,  TX  77001 

72.  Prof.  Steven  K.  Thompson . 

Statistics  Dept. 

Pennsylvania  State  Univ. 

326  Classroom  Bldg. 

University  Park,  PA  16802-2111 

. 1 

73.  Prof.  J.  W.  Tukey . 

1 

Statistics  Dept.,  Fine  Hall 

Princeton  University 

Princeton,  NJ  08540 

74.  Dr.  D.  Vere-Tones . 

1 

Dept,  of  Math 

Victoria  Univ.  of  Wellington 

P.O.  Box  196 

Wellington,  NEW  ZEALAND 

75.  Prof.  David  L.  Wallace . 

1 

Statistics  Dept. 

University  of  Chicago 

5734  S.  University  Ave. 

Chicago,  IL  60637 

76.  Dr.  Ed  Wegman . . . 

1 

George  Mason  University 

Fairfax,  VA  22030 

77.  Dr.  L.  Wein . 

1 

Operations  Research  Center,  Rm  E40-164 

Massachusetts  Institute  of  Technology 

Cambridge,  MA  02139 

78.  Dr.  Alan  Weiss . 

1 

Rm  2C-118 

AT&T  Bell  Laboratories 

600  Mountain  Avenue 

Murray  Hill,  NJ  07974-2040 
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79.  Prof.  Roy  Welsch . . . 1 

Sloan  School 

M.I.T. 

Cambridge,  MA  02139 

80.  Dr.  Raymond  S.H.  Yang . 1 


Colorado  State  University 

College  of  Veterinary  Medicine  and  Biomedical  Sciences 
Dept,  of  Environmental  Health 
Fort  Collins,  CO  80523 
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